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Abstract 

We study the effects of a smah density of holes, 5, on a square lattice antiferromagnet undergoing 
a continuous transition from a Neel state to a valence bond solid at a deconfined quantum critical 
point. We argue that at non-zero 5, it is likely that the critical point broadens into a non- Fermi 
liquid 'holon metal' phase with fractionalized excitations. The holon metal phase is flanked on 
both sides by Fermi liquid states with Fermi surfaces enclosing the usual Luttinger area. However 
the electronic quasiparticles carry distinct quantum numbers in the two Fermi liquid phases, and 
consequently the ratio Wms^QAp/S (where Ap is the area of a hole pocket) has a factor of 2 
discontinuity across the quantum critical point of the insulator. We demonstrate that the electronic 
spectrum at this transition is described by the 'boundary' critical theory of an impurity coupled 
to a 2+1 dimensional conformal field theory. We compute the finite temperature quantum-critical 
electronic spectra and show that they resemble "Fermi arc" spectra seen in recent photoemission 
experiments on the pseudogap phase of the cuprates. 
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I. INTRODUCTION 



There was a great deal of work on the dynamics of a single hole in a square lattice anti- 
ferro magnet, soon after the discovery of high temperature superconductivity in the cuprate 
compounds. It was demonstratedii2i'^i^»^i^'^>^ that a single hole moving in a Neel ground state 
has a finite quasiparticle residue, Z\ so a small density of holes, 5, are expected to form 
a Fermi liquid. This Fermi liquid state with Neel order will be the starting point of our 
analysis. Also, Shraiman and Siggia^*^ introducing a current- current coupling between the 
hole and the antiferromagnet which implied that a large spin S Neel state is unstable for 
certain parameter ranges to spiral spin ordering; we shall not be interested in this metallic 
spiral state here, although the Shraiman-Siggia coupling (in Eq. (12.121) below) will play a 
key role in our analysis. 

We begin with a 5* = 1/2 Neel state of an insulating antiferromagnet and imagine "turning 
up quantum fluctuations" by adding further neighbor or ring-exchange couplings so that 
there is a transition to a paramagnetic state in which spin rotation invariance is restored. 
Now add a small density of holes to this antiferromagnet. The main question we shall 
address is: what is the fate of the Fermi liquid Neel state across such a transition ? 
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Specifically, we consider the 'deconfined' quantum phase transition proposed in Ref. 
for an insulating = 1/2 square lattice antiferromagnetp^»i^'^»i^. This is a theory for a 
transition between a Neel state and a spin-gap state with valence bond solid (VBS) order 
(the latter state is spin rotation invariant, but breaks lattice symmetries by ordering of 
valence bonds). These two states break distinct symmetries of the Hamiltonian, and so 
cannot generically have a continuous transition between them in the Landau-Ginzburg- 
Wilson theory of phase transitions. However, such a transition is found in a 'deconfined' 
theory focusing not on order parameters but on fractionalized excitations and emergent 
gauge forces. The transition is tuned by the coupling s (which represents the strength of 
frustrating exchange interactions) — see Fig. [H Upon doping, we will argue that the most 
likely possibility is that the insulating deconfined critical point gets broadened into a novel 
non- Fermi liquid 'holon metal' phase, with no Fermi surface (shown shaded in Fig. [1]). 
The qualitative distinction between the Neel and VBS states survives also in the Fermi 
liquid states at 5 > (shown unshaded in Fig. [1]): we will show that a characteristic 
property (specified shortly) of the Fermi surface has a discontinuity in the limit 5 — >■ 
as s is scanned across the s = critical point of the insulator. We also compute finite 
temperature electronic spectra in the vicinity of the transition and find that they resemble 
"Fermi arc" spectra seen in recent photoemission experiments on the pseudogap phase of 
the cuprates^. 

There have been other discussions in the literature^ii^'22 of transitions between Fermi 
liquids, including proposals that there could be a continuous quantum transition with a 
discontinuous change in the shape of the Fermi surface (recent experiments^^ on CeRhlus 
are compatible with an abrupt or very rapid change in Fermi surface topology). This would 
require a sudden change in the Fermi surface as a function of s at a fixed non- zero value 
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FIG. 1: Schematic phase diagram. The elhpses represent spin singlet valence bonds. The coupling s 
tunes the insulator across the Neel-VBS transition, and 5 is the mobile hole density. The deconfined 
quantum critical point is at s = Sc in the insulator with 5 = 0. The vacancies ('holons') carry a 
gauge charge q = ±1 under an emergent U(l) gauge force. In the cartoons above, the reader can 
interpret g as a sublattice label. In the s < Sc, Neel phase, q determines the spin: a vacancy on an 
up (down) spin site carries net spin down (up) and so is equivalent to a charge e spin-1/2 hole. For 
s > Sc, the hole is a composite of a vacancy and a nearby unpaired spin with opposite q, moving 
by rearranging nearest-neighbor valence bonds; note that this motion preserves spin and sublattice 
quantum numbers separately (see also Ref. Il6l ). So there are twice as many states per momentum 
for a charge e spin 1/2 hole in the VBS state than there are in the Neel state. 



of S. We will argue that such a change is unlikely in our models, and the situation is as 
illustrated in Fig. [1], with an intermediate non- Fermi liquid phase. 

By an extension of arguments in early work— i^^i^i^^ it is expected that a significant 
portion of the phase diagram in Fig. [T] is unstable at low temperatures to superconductivity. 
We defer consideration of such superconducting states to future work, and limit ourselves 
here to the normal states. 

We now turn to a more detailed summary of our results. First, we discuss our results in the 
unshaded regions of Fig[Tl In these regions, we are adding a small density of mobile carriers 
to conventionally ordered insulators, and we obtain Fermi liquid phases with electron-like 
quasiparticles with a non- zero quasiparticle residue, i.e. Z ^ 0. Non- Fermi-liquid physics 
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FIG. 2: Momentum space Fermi surfaces in the Neel and VBS regions of Fig[TJ The fihed circles are 
the 4 Kp wavevectors, with Ki = {n/2a){l, 1), K2 = (vr/2a)(l, -1), K3 = -Ki, = -K2, with 
a the lattice spacing. The dashed line in the Neel phase indicates the boundary of the magnetic 
Brillouin zone. Only the Fermi surfaces within this zone contribute to the Luttinger counting, and 
so the area of each ellipse is Af = (27r)^(5/4. In the VBS phase, all 4 pockets are inequivalent, 
and so the area of each ellipse is Af = (27r)^(5/8. The dashed lines now show the reduction of 
the Brillouin zone due to the VBS order which appears at sufficiently low temperatures; "shadow" 
Fermi sufaces, with weak photoemission intensity (estimated in the text), will appear as reflections 
across these lines, and these Fermi surfaces are not shown. 



appears only in the shaded region. 

In the s < Sc Neel phase, we obtain a Fermi liquid statei'^i'^i^'^'^'^ with four hole pockets 
centered at the Kp = (7r/2a)(±l, ±1), where a is the lattice spacing, shown in Fig. [2l 
However, because of the halving of the Brillouin zone by magnetic order, only two of these 
pockets are distinct. After accounting for the two- fold spin degeneracy, we conclude that the 
area enclosed by each pocket is Af = (27r)^5/4. Another way of understanding this halving 
of the Brillouin zone (which is also indicated in the caption of Fig. [1], and discussed further 
in Section IIIip is as follows. We can consider the doped hole as a vacancy in the background 
of a Neel state. If this hole is to move without leaving a string of broken bonds^, it must 
preserve its sublattice label. However, because of the broken symmetry associated with the 
Neel order, the sublattice location is not independent of the spin of the vacancy, and two 
labels are really the same quantum number. 

Next, we discuss a small density of holes in the s > Sc VBS state. As we will demonstrate 
in Section HV] in this state the four hole pockets are no longer pinned at the Kp, but instead 
move a distance A away, as indicated in Fig. [2l This shift arises from the Shraiman-Siggia^ 
coupling. The value of A is determined by s — Sc, but is independent of S to lowest order in S. 
Consequently, for sufficiently small 6, the hole pockets do not intersect the reduced Brillouin 
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zone boundaries, associated with the appearance of VBS order. The four hole pockets 
therefore all contain distinct quasiparticles states, and after accounting for the two-fold spin 
degeneracy, we now conclude that the area enclosed by each pocket is Ap = (27r)^(5/8. As 
above, another interpretation of this result is indicated in Fig [1], and will be described in 
more detail in Section llVi the hole motion in the VBS state also preserves its sublattice 
index, but now the sublattice and spin labels are distinct quantum numbers. 

We can summarize the above statements into one of the main zero temperature results 
of this paper: 

Af 



Af 
lim — — 

<5^0 6 



2 X lim , 

<5-*0 6 



(1.1) 

S>Sc 

Note that on both sides of the equation, we are taking the limit 5 ^ at fixed s. Thus, 
although a characteristic feature of the Fermi surface (the ratio Af/S) changes discontin- 
uously as s crosses Sc, the Fermi surface itself is of vanishingly small size. The result in 
Eq. (11.11) does not constitute a discontinuous change in the Fermi surface in the sense of 
other proposals^ii^i^. 

Next, we address the evolution of the Fermi surface along the fixed 6 > line P in Fig[Tl 
Aspects of the physics will be addressed in Section |Vl and some questions are deferred to 
future work. The key issue is the fate of the monopoles in the U(l) gauge theory which 
describes the deconfined critical point of the insulator at s = Sc- The situation has been 
discussed at length elsewhere^ for 5 = 0: the monopoles are irrelevant at the s = Sc critical 
point, but are relevant for all s > Sc where they induce confinement and VBS order. There 
are two distinct interesting possibilities for the fate of this transition for 5 > (a first-order 
transitions is, of course, also possible): 

(i) The critical point at s = Sc, S = extends into a single line of deconfined critical points 
for 6 = 0, with Fermi-liquid physics on either side of the line. In this case, physics discussed 
above Eq. (II. ip will apply also for 6 > 0, and there will be a discontinuous Fermi surface 
change across a continuous transition, realizing scenarios postulated in previous work— i^"^. 
However, we will argue in Section |V] that this possibility is unlikely in the present model. 
(a) The deconfined critical point at s = Sc, 5 = broadens into a deconfined phase, 
extending over a finite range of s values (shown as the shaded region in Fig. [1]), with novel 
non-Fermi-liquid physics. The monopoles, which were suppressed by the gapless critical 
modes at the single point s = Sc in the 6 = insulator, are now suppressed^^ over a finite 
range of values of s by the additional gapless excitations associated with mobile charge 
carriers for 6 > 0. Arguments supporting this possibility will appear in Section |Vl For 
reasons which will become clear in Section |Vl we will refer to this phase as a 'holon metal'. 
This holon metal phase has no electron-like Fermi surfaces, and hence Z = 0. 

It is convenient to think about the phase diagram in the plane where /i is the chemical 
potential. For a range of fi in the Mott insulator the ground state does not change for fixed 
s. The point s = Sc in this Mott insulator is the deconfined quantum critical point and 
its low energy theory is a conformal field theory (CFT). Below some critical /ic the hole 
density 6 holes will increase from zero as some power of /i — fic- It should be clear that 
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the critical point at s = Sc, /U = /Xc is an unstable multicritical point which controls the 
properties of the holon metal phase at low doping. In this perspective^, we view /ic — yU as a 
relevant perturbation to the field theory describing this multicritical point. So at an energy, 
temperatures, or wavevectors higher than a low energy scale determined by 5, the physics is 
described by the response of the CFT of the deconfined critical point of the Mott insulator 
to an infinitesimal hole density. We will argue in Section |V] that this response is associated 
with a 'boundary CFT' describing a quantum impurity in the bulk, 2+1 dimensional CFT. 
The result in Eq. (11. ip describing the Fermi liquid phases outside the shaded region of Fig. [1] 
can also be viewed as a consequence of the structure of this boundary CFT. 

Section IV Al will explore practical consequences of the boundary CFT by presenting nu- 
merical results on the form of the finite temperature electronic spectral weight. This will be 
carried out using the full action presented in Section [Til including terms that are formally 
corrections to scaling to the boundary CFT, but are important for experimental compar- 
isons. We will find that the results resemble recent photoemission measurements on the 
pseudogap phase of the cuprates^; in particular, the length of the "Fermi arcs" is roughly 
proportional to the temperature. Furthermore, we find that the Shraiman-Siggia^ term, 
which was responsible for the shift in the centers of the hole pockets away from the Kp in 
the zero temperature VBS state (Section [IV] and Fig. [2]), also shifts the centers of the finite 
temperature "Fermi arcs" from the Kp, as is seen experimentally. This feature is not present 
in a recent theory based upon a ' staggered- fiux' spin liquid^. 

We will begin in Section [III by presenting the field theory which describes the CFT at 
s = Sc, S = and the perturbations which induce mobile carriers for 6 > 0. The Fermi liquid 
phases of this field theory will be discussed in Sections IIIII and IIVI for the Neel and VBS 
states respectively. The non-Fermi-liquid holon metal phase, and the finite temperature 
quantum criticality is described in Section [Vl and we conclude in Section [Vll 

II. FIELD THEORY 

There is a vast literature on the theory of the lightly doped antiferromagnet. Two main 
approaches have been taken. The 'spin-fermion' models^, which are extrapolations of the 
weak-coupling Hartree-Fock theory of the Hubbard model, offer a convenient description of 
the Fermi liquid states with Neel order. However, these models cannot describe the transition 
to the VBS state, and so are not suitable for our purposes. The second approach,-!^ departs 
from the spin- wave theory of the insulating antiferromagnet, and describes the coupling of 
the holes to a spin deformations of the antiferromagnet. In modern terms, this is formulated 
as 'chiral perturbation theory' in which the broken symmetry of the Neel state allows one 
to place constraints on the couplings between the mobile holes and the 'chiral' spin wave 
excitations^^. However, the non-linear realization of the SU(2) spin rotation symmetry in 
the Neel state makes this approach cumbersome. For our purposes, a central shortcoming 
is that this approach is not able to describe the fermion spectrum in a state in which SU(2) 
symmetry is restored and the magnetic Brillouin zone expands to the full Brillouin zone of 



6 



the square lattice: Bloch reflections across the magnetic Brillouin zone boundary do not 
disappear at any order in chiral perturbation theory. It is essential to considerations of 
Fermi surface topology that such effects are properly accounted for. 

We will instead use a method in which the spin SU(2) symmetry is linearly realized 
and is kept explicit at all stages. We also want to explicitly preserve the full space group 
symmetry of the square lattice. In the insulator, this requires starting from a 'spin liquid' 
state. Although there is no gapped spin liquid insulator in FiglH the VBS state descendsiii^ 
from an instability of a particular U(l) spin liquid,— and the latter will serve as our starting 
point. Also, the CFT at s = Sc, S = 0, defines an 'algebraic spin liquid', and this is obtained 
as a gapless critical point in the theory of the same U(l) spin liquid. 

It has been emphasized by Wen^ that an essential characteristic of any spin liquid is its 
PSG: the projective symmetries of the various fields under elements of the square lattice 
space group. Because the fields carry charges under a gauge group which characterizes the 
spin liquid, they need not be invariant under such transformations — they are only required 
to be invariant up to a gauge transformation, and this is sufficient to preserve square lattice 
symmetry in all observables. Here we will show how the PSG can be extended from the 
insulator to the doped state, where it also places crucial constraints on the low energy 



effective theory. Another PSG-based analysis of a doped Mott insulator is in Ref. |33| for 
case where the spin liquid has fermionic spinous. Here we focus on the vicinity of the 
Neel-VBS transition where, as we review below, the spinous are bosonic. 

The motion of charge carriers doped into an insulating square-lattice quantum antiferro- 
magnet is conventionally described by the "t-J" model, 

Ht-j = - "^Ujicfcja + h.c.) + ^ JijSi ■ Sj + . .. , (2.1) 

where creates an electron with spin a on the sites i of a square lattice and Si = 
h^ap^'i'^ ^a^il^^ with c? the Pauli matrices. In addition, the constraint ^Q,c°^CjQ, < 1 is 
enforced on each site, modeling the large local repulsion between the electrons. It is im- 
portant to note that our results are more general than a particular t-J model, and follow 
almost completely from symmetry considerations. The ellipses in Eq. (12. ip additional short- 
range couplings which preserve square lattice symmetry and spin rotation invariance. The 
coupling s is some axis in this multidimensional parameter space which accesses both the 
Neel and spin liquid phases in the insulator. 

We will analyze Ht-j by choosing a representation of the electron operator which obeys 
the constraint on each site, and most conveniently accesses the U(l) spin liquid describing the 
Neel-VBS transition. This is achieved by introducing neutral Schwinger bosons to represent 
the spin degrees of freedom, and spinless fermionic holons which track the charge. Although 
this is not essential, we will use distinct representations on the two square sublattices because 
this simphfies the mean-field structure of the spin liquid solution. On one sublattice of the 
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square lattice we write the electron operator, 

Ca = bjl (2.2) 

where ba are canonical Schwinger bosons and /+ are canonical fermionic holons. The con- 
straint on each site now becomes 

+ = (2.3) 
On the other sublattice, we use bosons which transform as a conjugate respresentation 

Ca = Bafsffl (2.4) 

with a similar constraint. We define the antisymmetric tensor by e^i = 1 and Ea/s = s°'^ ■ 
Notice that the representations in Eq. (12.2112. 41) have a U(l) gauge redundancy. The structure 
of the mean-field theory described below instructs us to take the continuum limit so that 
the U(l) gauge charges of ha and /+ are +1, while those of the 5" and /_ are —1. 

First, we recall the mean field theory^ of Ht_j in the insulator in the above represen- 
tation. Here we can ignore the holons f±. The exchange interactions are quartic in the 
Schwinger boson operators, and we decouple these by a Hubbard-Stratanovich transforma- 
tion. This yields the mean field Hamiltonian for the U(l) spin liquid 

Hj = -QY, ^-^J + H-c. + A 5^ 6f 6,, + A ^ 6] j;, (2.5) 

{ij) « i 

where i is restricted to be on one sublattice, and j on the other, and Q, A are positive 
constants to be determined by solving the mean-field equations. Upon adding frustrating 
exchange couplings to ifj, the mean field parameters will vary, allowing us to access both 
the Neel and spin liquid phases. For sufficiently large A, the boson spectrum is gapped, and 
this describes the spin liquid (s > Sc). As A is decreased, the bosons eventually condense, 
leading to a transition to the Neel state. We now want to obtain a continuum field-theoretic 
representation of this quantum phase transition of the insulator. This can done by carefully 
taking the continuum limit of the low energy excitations of Hj, and then including fiuctua- 



tions about the mean-field saddle point: such a procedure is described in detail in Ref. |30 . 
Here, we show that the same answer can be rapidly obtained by a PSG analysis. For this, 
we need to understand how Hj preserves the full square lattice symmetry even though it 
treats the two sublattices inequivalently. So we consider the transformations of the opera- 
tors under all the space group operations of the square lattice, and also under time-reversal. 
These are listed in Table [H From the PSG of the Schwinger bosons, and using Eqs. (12.2112.41) . 
we can immediately deduce the PSG of the lattice fermion fields f±. These are also shown 
in Table [11 Note especially the sign in the last row, which is a consequence of e^'^Sp^ = — 

We now use the PSG to obtain the continuum theory of the Neel-VBS transition in the 
insulator. The boson spectrum of Hj shows that the minimum energy excitations are near 
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TABLE I: PSG transformations of the lattice fields under square lattice symmetry operations. 
Tx'. translation by one lattice spacing along the x direction; R^^'- 90° rotation about a dual 
lattice site on the plaquette center {x ^ y,y ^ —x); I^^^^'- reflection about the dual lattice y axis 
(x — > —x,y y); T: time-reversal, defined as a symmetry (similar to parity) of the imaginary 
time path integral. Note that such a T operation is not anti-linear. The transformations of the 
Hermitian conjugates are the conjugates of the above, except for time-reversal of fermions^^. For 
the latter, f± and are treated as independent Grassman numbers and T : — > —f±. 

zero momentum, and so we may take its low energy limit by a naive gradient expansion. 

— cx. 

We can proceed in terms of the 6^ and the h , but it is convenient to introduce the linear 
combination 

~ + hi (2.6) 

because the orthogonal linear combination can be integrated out--. The PSG of the 'spinous' 
Za follows from Table [T] and is listed in Table [Tll Also important for the low energy theory is 
a U(l) gauge field, (/i = r, x,?/ is a spacetime index), which represents the fluctuations 
of the phase of Q and of A. The Za carries charge -|-1 under the gauge interaction. The 
explicit PSG of the can be found elsewhere^, and we do not list it explicitly here because 
gauge-invariance is sufficient to determine the terms. Writing down the most general 
action invariant under the PSG we obtain the theory of the low energy excitations of the 
spin liquid: / (PrdrCj with 

A = |(9^-^A^)^j2 + s|z,p + |(|z„P)' (2.7) 

Here s is the tuning parameter in Fig. [H we have set a spinon velocity to unity, and u is 
a spinon self-interaction. For s > Sc, this field theory is in a spin SU(2)-invariant U(l) 
spin liquid phase, with spinon excitations above an energy gap A. This gap vanishes as we 
approach the critical point as A ~ (s — ScY, where u is the correlation length exponent 
of the CFT describing the deconfined critical point. For s < Sc, the Za condense into a 
Higgs phase, which is the Neel state of the antiferromagnet. The Neel order parameter is 
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TABLE II: PSG transformations of the fields in the continuum theory under square lattice sym- 
metry operations. 

N = z°'*a^Zi3, and this order vanishes as {sc — 

The term Cm describes the dynamics of monopoles in the compact U(l) gauge field A^. 
This term has been discussed at length elsewhere,—'^ and we will not reiterate the details 
here. For our purposes, the important conclusions of these earlier analyses are: (z) monopoles 
are suppressed for s < Sc, in the Neel phase and at the deconfined quantum critical point; 
(m) monopoles proliferate in the s > Sc spin liquid phase, inducing confinement of spinous 
and the appearance of VBS order. A key fact^'^ which leads to the latter conclusion is 
that Berry phases of the underlying antiferromagnet endow the monopoles with a non-trivial 
PSG^^. Indeed, this PSG allows us to identify^ the monopole creation operator with the 
VBS order parameter, $vbs- This order measures modulations in the density of singlet 
valence bonds on the square lattice, and is defined by 

$VBs(r1 = {~lf^Sr-S,+,+l{-lY-ySr-Sr-+y, (2.8) 

where x, y are the unit lattice vectors. The proliferation of monopoles implies that for s > Sc, 
('^'vBs) 7^ 0. This order vanishes as we approach the critical point, ($vbs) ~ (s — Sc)'^'^^^, 
and the exponent is determined by the scaling dimension of the monopole creation operator. 
The representation in Eq. (12. 8p also allows us to determine the PSG of $vbS) and hence of 
the monopoles, and this is also listed in Table [Tll 

Now let us move to the doped antiferromagnet, and examine the nature of the low energy / 
motion. The primary input we need from the microscopic physics is the form of the dispersion 
of a single / fermion in the Brillouin zone in the background of the U(l) spin liquid just 
described. For t ^ J, this leads to a strongly-coupled problem and accurate analytic results 
are not possible. However, numerical data is available on the dispersion spectrum in the 
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Neel phasei'^i'^'^i^'^, and we know that for a model with nearest-neighbor hopping, t, a single 
hole has its dispersion minimum at the four Kp points shown in Fig. [2l We will proceed 
here using the results of a self-consistent analysis of a single / excitation above the U(l) 
spin liquid state as carried out in previous work-i^. This mean field analysis is reviewed in 
Appendix |A] in our notation. As we will see in Section IIIIl the resulting theory derived here 
predicts a hole dispersion in the Neel phase consistent with the numerical results^'^'^'^'^i^. We 
also note that for suitable second-neighbor hopping, t', the hole dispersion in the Neel phase 
shifts its minima away from the Kp points to (7r/a)(0, ±1), (7r/a)(±l, 0): this is believed 
to be relevant for the electron doped cuprates^^. The PSG analysis presented below will 
require modifications of this case, but will not be discussed here. 

We already know that the fermionic holons, /, carry a U(l) gauge charge g = ±1. The 
analysis in Appendix lAl shows that they also acquire an additional 'valley' quantum number 
because their dispersion minimum is not at the center of the Brillouin zone; we will keep 
track of this at the cost of some additional bookkeeping. There are 2 distinct valleys, with 
labels V = 1,2, and we choose a gauge in which their minima are at wavevectors Ki^2 shown 
in Fig. [21 Summarizing, we have 4 species of spinless fermions fqy, with U(l) charge q = ±1, 
and valleys v = 1,2. The PSG of the /g„ can now be deduced from the PSG listed in 
Table [H We need to modify the transformations here by additional factors arising from a 
factor e*'^"'^, associated with shift of the dispersion minimum to a finite wavevector. Such a 
procedure leads the results in Table HTl 

Before we turn to the effective action for the it is useful to write down an expression 
for the physical electron operator. This will allow to compute observable properties of our 
theory, such as the location of the Fermi surfaces. This expression can be obtained by 
requiring a spinon-holon composite to be gauge invariant, and to transform under the PSG 
like a physical electron. Because the fgy fermions reside near the Kp momenta defined in 
Fig [21 while the spinous are near zero momentum, it is convenient to decompose the charge 
— e, spin-1/2 electron annihilation operator Ca{r) into electron-like components, "^pa, near 
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^e*^-%,(f). (2.9) 

p=i 



Now the PSG requirements yield unique bilinear combinations for the 
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^2,4a = zJl^±e^pz''*fl^. (2.10) 

We finally turn to the effective action for the fq^ holons. First, the terms bilinear in the 
/, invariant under the PSG, are: 

^f = Y.fl - iqAr - /z - ^^^^^) (2.11) 
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where j extends over x,y, mix = fT^2y and m2x = ^ly are the mass of the elhptical hole 
pockets, is the hole chemical potential, and the x and y directions are rotated by 45° from 
the principle square axes. The PSG prohibits a linear derivative term in Eq. (12.111) . and so 
the dispersion is an extremum at zero momentum; this pins the Fermi surfaces to be 
centered at Kp in the Neel phase (Fig. [2]), but not, as we will see, in the VBS phase. 

Next, we include coupling between the holons and spinous invariant under all PSG oper- 
ations. There is an unimportant scalar coupling Yliaqv \'^af' fqvfqv) but the more important 
term is 

C-c = «Ae"^ \^fllf-lZadxZf3 + fl2f-2ZadyZp^ + C.C., (2.12) 

equivalent to the dipole coupling introduced by Shraiman and Siggia^, arising from the 
hopping of electrons between nearest-neighbor sites. 

Finally, we have to consider couplings between the holons and the U(l) gauge dynamics. 
The minimal coupling to the gauge field is already included in Additionally, we have 
to allow for couplings between the monopoles and the holons. These are also specified by 
the PSG. First, we search for terms which involve $vbs and a bilinear of the fgy fermions. 
However a detailed analysis shows that there are no such terms which are invariant under 
all the PSG operations; this is seen by first listing the bilinear invariants under I!^'^^^ 
(under which $vbs is invariant), and then noting that their transformations under are 
incompatible with those of $vbs- Consequently, the coupling between the VBS order and 
the fermions will be weaker than might have been initially expected, and will vanish faster 
than ($vBs) ~ A^^bs/i' as A ^ 0. The simplest non-vanishing coupling turns out to require 
the full electron operators in Eq. fl2.10p . This has the form 

-^VBS = A VBS VBS 



+ (^1^2 - ^2^1 + ^1^3 - ^3^4) +C.C. (2.13) 

Our complete field theory of the doped antiferromagnet is now specified by the Lagrangian 

C = Cj + Cf + C, + CyBs- (2.14) 

The following sections will describe its properties in the various states in the limit of small 
hole density 6. 



III. DOPING THE NEEL STATE 

As noted earlier, there is much earlier work on the Fermi liquid state obtained by doping 
a Neel state^'^'^'^'^'^'^. Here we present an instructive argument which shows how these 
results are reproduced in the present field-theoretic context. In particular, we will show that 
the gauge charge index q and the spin quantum number are tied to each other in the Neel 
phase, and also obtain the results quoted earher on the volume of the hole pockets. 
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Consider the s < Sc Neel phase of Cz with (za) ^ 0. By the Higgs mechanism, the 
condensate of renders the Ay, massive, and so we can safely integrate the out. A 
crucial point is that the Higgs mechanism also ties the U(l) gauge charge q of the holons 
to the spin quantum number 5*2 along the Neel order (because of the broken symmetry, Sx 
and Sy are not conserved); more precisely Sz = q/^, and so the fqy carry all the quantum 
numbers of the electron. By spin rotation invariance, we can always rotate the Za condensate 
(and without any rotation of the spinless fqy) to produce a Neel order iV" = z°'*a'^^zp (where 
(7° are the Pauli matrices) polarized along the (0, 0, 1) direction. Now examine the response 
of the theory to a uniform magnetic field H applied along the z direction. Under such a 
field, the only change in the action is that^ 

\{dr - iAr)za\^ ^ {{dr - lA, - Ha'/2)z^) X {{dr + lAr + Ha'/2)zl) . (3.1) 

Choosing {za) = ^J\s — Sc\/u5a^, we obtain a term in the Lagrangian 

£ = • ■ ■ (|s - s,\lu) [lA, - (i//2) V]' + • • • (3.2) 

which gaps the A^ photon. Note, however, that Aj. fluctuates about a non-zero value that 
is set by H. Integrating out the A^ and then evaluating = 6S/6H\jj^q, we obtain an 
expression for the magnetization density Mz 

Mz = \Y.^fU^ + --- (3.3) 

qv 

where the ellipses represent an additional term which measures the magnetization of the 
spin waves. This establishes, as claimed in Fig. [H that Sz = q/2 for the fermions in the Neel 
phase. 

Now add a single hole and examine its coupling to the spin excitations of the Neel state. 
The Afj, fluctuations (which also represent spin waves in the Neel state) are expected to 
strongly renormalize the effective mass of the hole. The Cc term in Eq. (12.121) does not 
modify the fermion dispersion at first order in A, but does lead to a finite correction to the 
fermion mass at second order. Consequently fqy fermions have a stable dispersion minimum 
at zero momentum; or in other words, by Eq. (I2.10p . the dispersion of the physical electrons 
has minima which remain pinned at the Kp. We will see that the influence of the Cc term 
is quite distinct in the VBS state. 

Now we consider the Fermi liquid state obtained with a total density 6 of the four fermion 
species fqy. Each will form a separate Fermi surface containing 5/4 states. From Eq. ( ]2.10p 
we deduce that there are 4 hole pockets centered at the Kp wavevectors, each enclosing the 
area Ap = (2vr)^((5/4), as shown in Fig. [2l The caption of Fig. [2] shows that the same area 
is obtained in direct Fermi liquid counting of electrons within the magnetic Brillouin zone. 
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IV. DOPING THE VBS STATE 



This section will establish one of the new results of the paper: the nature of the Fermi 
liquid state obtained by doping the VBS state. We will show that the electronic quasipar- 
ticles have a dispersion which is not centered in the Kp, and consequently the volume per 
hole pocket chang shown in Fig. [2J 

First, we recall some essential characteristics^!^ of the U(l) spin liquid for s > Sc, 
S = 0. This state has a "photon", representing a single excitation associated with 
rearrangements of valence bonds. The effective action of this photon can be estimated in 
the large N limit, where the spin index a = 1 . . . N. In this limit we integrate out the Za 
quanta in the action Cj, and obtain an effective for the A^. Details appear in Appendix [Bl 
and at scales larger than the spin correlation length, ~ 1/A, we have the Lagrangian 

= ^ ie,,xd.A,f (4.1) 

As shown in Appendix IB 11 in the non-relativistic limit of slowly moving excitations, this 
photon mediates a logarithmic "Coulomb" interaction V{r) between charged particles^^i^: 

V{r) = ^^\n{rA) , (4.2) 

which also appears in Eq. ( 1B11I) . The divergence of V{r) as r ^ oo implies that only 
configurations which are net charge zero have finite energy. Of course, at sufficiently large 
r, we also have to consider the modifications due to monopolesii*^^^ which change the 
mediated interaction from logarithmic to a linearly confining one. However, this does not 



happen until a confinement length scale denoted ^vbs in Ref. [Ill, whose divergence as A ^ 
is determined by the scaling dimension of the four-monopole operator, so that ^vbs ^ 1/A. 
So there is substantial scale over which monopole effects can be neglected, and we can work 
with the V{r) in Eq. flBTTl) . 

Now let us examine the spectrum of a single / hole in the VBS state. We will see that its 
key quantum numbers are determined on a scale smaller than ^vbs^ and so we can work with 
the non-compact gauge field obeying the action in Eq. (14.11) . As shown in Appendix IB 11 
the coupling to A^- causes this holon to have an 'electrostatic' self energy which diverges 
logarithmically with system size^, and so a spin-singlet single holon state is not stable. 
Rather, the holon will peel off a single spinon from above spinon gap, and form a S" = 1/2, 
charge e bound stat^. This bound state is electron-like and is neutral under the U(l) 
gauge force. A finite density of such bound states can then form a Fermi surface with charge 
e, S" = 1/2 quasiparticles. We also have to consider the pairing between holons with opposite 
U(l) gauge charges, induced by the gauge force (this attractive force must be balanced 
against the repulsive Coulomb force associated with the physical electromagnetic charge e 
of each holon). This is expected to lead to superconductivity at low enough temperatures, 
but we will not discuss this here; our focus is on the normal state. 
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We now discuss the wavef unction and quantum numbers of the holon-spinon bound state. 
This can be addressed in a non-relativistic theory of holons and spinous interacting via 
a Coulomb force, as discussed in Appendix IB 1[ The force in Eq. fIBlip will bind the 
spinons/antispinons with the holons into gauge neutral combinations. For each holon species 
fqy, there are two such bound states, because the spinon can have spin up or down. So there 
are a total of 8 electron-like bound states: this is the crucial difference from the electron- 
like bound states in the Neel state, where there were only 4 electron-like quasiparticles. 
Another key difference is that the Shraiman-Siggia term, Cc, mixes these 8 bound states, 
with significant consequences for the structure of the Fermi surfaces. 

The form of Cc in Eq. f l2.12p makes it clear that the bound state of a /+i holon with a 
spinon will be mixed with the bound state of a f^i holon with a spinon (parallel considera- 
tions apply to the f±2 holons). Before writing the Schrodinger equation for this bound state, 
we need to decompose the relativistic field Za into non-relativistic, canonical boson fields Pa 
(the spinon) and h"' (the anti-spinon). At low momenta, this is done by the parameterization 
(as in Eq. M ) 

{Pa + Bo^ph"^) (4.3) 



V2A 

We now consider the two-component wavefunction describing the Cc induced mixing 
between the bound state of /+i and ha and the bound state of /_i and Pa- Let us label 
the co-ordinates of the holon by fh and that of the spinon by r,,. Then the bound state 
wavefunction has a doublet structure between these states 

*(f„?.) = e-«(^+«) (4.4) 
and is expressed in terms of the center-of-mass and relative co-ordinates 

^1 — A 

f = fh-fs, (4.5) 

where A is the mass of the spinon in the non-relativistic limit, and the holon masses 

in Cf. Similarly, we can also define the total masses and the reduced masses 

Mj = m„j + A 
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For completeness, we also note the operator representation of this state: 



\aP+) 
\aP-) 



dr-hdrhe'^'^<P+ (f ) /j^ (f,) /i"t (^-;) | g) 
rfr-i(irle^^V(r)/!i(r,K^(f,)|0) 



(4.7) 



Using a reduction closely related to that presented in Appendix IB 1\ we deduce that the 
Schrodinger equation obeyed by 0±(r) is then 



P-2 
3_ 

2Mt 



+ V{r) 






(4.8) 



We first examine this equation at A = 0. For small spin gap A, the U(l) gauge potential 
potential in Eq. fIBlip . V{r) = {12A/N) In(rA), binds the holons and spinous over a length 
scale ~ 1/A, with 0±(O) = 0(0) ~ A. It is also clear that P dependence of Ep is only in 
the center-of-mass kinetic energy. The specific information we need from Eq. (14. 8 p is the 
modification of the dispersion Ep induced by Cc- This can be estimated in perturbation 
theory in A using the matrix element {aP + \Cc\c(P—) = —X\(f){0)\'^Px/Mx (we used the fact 
that 9^.0(0) = 0), and leads to 



E, 



Eo + 



P-2 

J 

2M-7 



± 



AR. 



(4.9) 



with eigenmodes which correspond precisely to the electron states in Eq. (I2.10p . So the 
bound state dispersion has a minimum at = ±A where A = A|0(O)p ~ A^. This is 
responsible for the shift in the centers of the elliptical hole pockets shown in Fig 2. 

The usual counting argument now allows us to deduce the area of each hole pocket. There 
are 4 inequivalent pockets, and a factor of 2 degeneracy for spin; so Af = (27r)^((5/8). The 
factor of 2 difference from the result obtained in the Neel phase is one of our key results. 

The state we have described so far is actually not a conventional Fermi liquid: the area 
enclosed by its Fermi surface is not the same as that of a non-interacting electron gas at 
the same filling and with the same size of unit cell. In the non-interacting case one would 
have a Fermi surface area (27r)^(l — 6)/2. In our case, where the holes were doped on a 
background spin liquid state with a gapless photon excitation, we obtain that the electron- 
like area enclosed by the Fermi surface is (27r)^(l — 6/2) . This state is a fractionalized 
Fermi liquidr^, obtained here in a single band model, in contrast to its previous appearance 
in Kondo lattice models. 

The instability of the spin liquid to confinement and VBS order induced by monopolesi^i^, 
and the associated halving of the Brillouin zone will finally transform the state into one 
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obeying the conventional Luttinger theorem, with an electron-hke area enclosed by the 
Fermi surface of (27r)^(l — 5)/2. There will be Bragg reflection across the reduced Brillouin 
zone boundaries. However, because of the shift in the minimum of the holon-spinon bound 
state dispersion discussed above, this mixing is negligible as long as < A ~ A^, and so 
can always be neglected as 5 — 0. This establishes the Fermi surface structure claimed in 
Fig.H 

At larger 6, the hole pockets will cross the dashed lines in Fig. [21 and Bragg reflections 
will split the Fermi surfaces. We can analyze this magnitude of the mixing matrix element 
in the Hamiltonian, using £vbs in Eq. fl2.13p . Using Eqs. fl4.3p and fl4.4l) . we obtain a matrix 
element of order 

AVBS($VBS)^^ ~ AvBsA^+^VBs/^. (4.10) 

As expected, this does vanish faster than A^^^^/'', and is responsible for the weak Bragg 
reflection across the reduced Brillouin zone boundaries of the VBS state in Fig. [2l The 
splitting of the Fermi surface is negligible provided the matrix element is smaller than the 
hole kinetic energy, or 6 > A^+^^bs/;^. Over such a possible regime of larger 6, the basic 
pictures of Figs. [T]and[2]continue to hold. The photoemission intensity of the "shadow" Fermi 
surfaces noted in Fig. [2]is proportional to the square of the matrix element or ~ /\2{i+/3vbs/'^)_ 



V. HOLON METAL AND QUANTUM CRITICALITY 

To complete the picture, let us address the physics of the shaded region in Fig. [H Ne- 
glecting the gauge fluctuations, the holons and spinous are independent. The four holon 
species then form their own Fermi-surfaces centered around -fi'i, 2,3,4- In the Neel phase, the 
spinous condense leading to long-range order. In this case the physical electron fields are the 
holons On increased doping, Neel order is destroyed, which corresponds to the development 
of a gap for the spinous. Now we have a finite density of holons that we expect will form 
a 'holon metal' with a holon Fermi surface. Including gauge fluctuations we have a theory 
of a finite density of holons fg^ interacting with spinous via a U(l) gauge force A^. At 
sufficiently long scales, the holons will 'Thomas- Fermi' screen that longitudinal forced, 
and so obviate the binding into gauge neutral combinations. Consequently the spinous and 
holons remain as relatively well-defined excitations, and we expect to enter a fractional- 
ized holon metal phase. This screening can be prevented if the spacing between the holons 
(~ 1/v^) is larger than the holon-spin binding length ~ 1/A, or 6 < A^; this fixes the 
boundary of the shaded region in Fig. [H Alternatively, the same results follows from our 
argument in Section [T] that the scaling dimension of 6 is 2. The shaded region will exclude 
the unshaded region of the VBS phase where the Fermi surface splitting is negligible (dis- 
cussed in the previous paragraph) provided /?vbs > ^, an inequality that holds at least for 
large N. Note that it is the competition between the holon kinetic energy and the gauge 
mediated attractive interaction with the spinous that determines the relative stability of the 
holon metal and the VBS phases. 
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It is this fractionalized phase that we propose as a candidate for the pseudo-gap phase 
of the cuprates. Due to the presence of the gapless gauge field the low energy properties of 
the holon metal will be modified in well-known ways from that of a non-interacting holon 
picture, and we will not repeat these further here. We emphasize again here that the "holon- 
metal" phase should be stable to confinement due to monopoles following previous screening 
arguments-^£. The metallic holon-metal phase will however inevitably be unstable at low-T 
to pairing and hence superconductivity. 

We now describe the criticality of the hole spectrum at s = Sc- We assume here an 
observation scale (frequency (u) or temperature (T)) large enough, or a 5 is small enough, 
so that the holes can be considered one at a time. A key observation about a single hole is 
that its quadratic dispersion {i.e. the terms proportional to l/{2myx,y) in ^Cf) is an irrelevant 
perturbation on the quantum critical point of which involves excitations which disperse 



linearly with momentum; a similar observation was made in Ref. |40| for Landau-Ginzburg- 
Wilson critical point of 0(3) model. Consequently, the hole may be considered localized, 
and its physics is closely related to the single impurity in a spin liquid problem analyzed 



in Ref. |37|. Here we are interested in the single hole Green's function and this requires 
the overlap between states of the spin liquid with and without the impurity. This can be 
computed by analyzing the quantum critical theory of coupled to one holon localized at 
r = 0, represented by the r-independent Grassmanian /(r): 

Sg, = j d^rdrC, + J drf^ + - tAr{r = 0, r)^ /. 

Here eo is an arbitrary energy fixing the bottom of the holon band, and, following earlier 
arguments^I, the only relevant coupling between the 'impurity' holon degree of freedom and 
the bulk degrees of freedom of is the gauge coupling associated with A^-. The spectral 
function of a physical charge e, S = 1/2 hole is given by the two-point correlation of the 
composite operator -Zq/^- If this operator has scaling dimension r]h/2, then the universal 
critical hole Green's function, Gh, is independent of wavevector and obeys the scaling form 

= T-^^-^'^^^{h{u - eo)/kBT) (5.1) 



where $ is a universal scaling function. At T = 0, we obtain an incoherent spectrum 
associated with the power-law singularity Gh ^ {uj — eo)~^'^^'\ We have computed rjh by 
a standard 1/A^ expansion for Gh under the action Sgc, the computation is described in 
Appendix IB 2\ and the result is 

Another perspective on the exponent rjh is obtained by mapping it to an observable in 
the lattice non-compact CP^ model studied by Motrunich and Vishwanath^^. On a three- 
dimensional cubic lattice with spacetime points j, the model has the complex spinor fields zja 
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on each lattice site, and a vector potential Ajfj, on each link extending along the /i direction 
from site j. The hole Green's function at imaginary time Mr (in units of the lattice spacing) 
is then given by the two-point correlator of the Zja in the time direction along with an 
intermediate 'Wilson line' operator (representing the contribution of the / fermion) which 



renders the correlator gauge invariant (see also Refs. I4llj42l ): 



GhiMr) = Zj* exp I 2 ^ Aj+nf,T I Zj+MrT,a 

' NCCP 

this is expected to decay as e~'^°^^ /M^'^ for large Mr at the quantum critical point. The 
continuum limit of the above correlator was computed by Kleinert and Schakel^ in the 1/N 
expansion: their result agrees with Eq. (15. 2p . 



A. Photoemission spectra 

The result Eq. (15. ip describes universal component of the electronic spectrum in the limit 
where "irrelevant" terms such as fermion dispersion and Cc in Eq. (12.121) are neglected. With 
an eye towards comparison with experiments!^, this subsection will present numerical results 
in the full Brillouin zone for a range of T with these effects taken into account. 

Formally similar computations have been carried out by Weng et al.^'^ at T = 0. However, 
they applied their result to the Neel phase, whereas we have argued that the Neel phase 
has a conventional Fermi liquid character. They also do not obtain the shift of the spectral 
weight maxima away from the -ft'i, 2,3,4 points in the Brillouin zone. 

We are interested in the electron spectral function that would be measured in a photo- 
emission experiment. The measured photo-emission intensity is related to the spectral den- 
sity of the physical electron Greens' function. 

G^^(r,r) = -(r[c.4r)4(0)]) (5.3) 

In order to re-write this expression in terms of the z and / fields we need an expression for 
the physical electrons in terms of these fields. As shown earlier, we can derive such a relation 
from the projective symmetry group (PSG) transformation properties of the various fields 
Eq. (ElOD, 

Cra = 6 ^ [Zaf\.i ~\~ ^afB^^ /— l] 

+ e''''-'-[zjl^-e^pz^*fl^] (5.4) 
Using this expression, we can calculate the contributions to Eq. (15. 3p under the C = Cz + 
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(A) 0{X^) 



(B) 0{X') 




FIG. 3: Diagrams contributing to A^{k,uj) at (A) zeroth order and (B) first order in Eq. (j2.12p . 
The wavy (dashed) hnes are the spinon (holon) propagators. The shaded vertex corresponds to A. 

In a "cartoon" mean-field state for the holon-metal phase we neglect the gauge fluctu- 
ations. The holons and spinous are independent (we temporarily also set A = 0, but shall 
calculate the effect of its presence shortly). As described above in the Neel state the physical 
electron fields are the holons and thus in this phase the the photo-emission response should 
have coherent peaks exactly where the holons do, i.e. (±7r/2, ±7r/2). This is consistent with 
photo-emission experiments in the doped cuprates with long-range Neel order present^. On 
increased doping and entering the 'holon-metal phase with a finite spin gap, the spinous are 
not condensed and the associated spinonic fluctuations play a crucial role in the measured 
electronic spectra. We will calculate the spectra within this mean field picture and simply 
convolve the free spinon and holon Greens functions to get the electron Greens function. 
Qualitatively the neglected gauge fluctuations are expected to have two effects. First the 
scattering of the gapless holons with the gapless transverse gauge fluctuations leads to a 
singular frequency dependent self-energy for the holon Greens function. Upon convolving 
with the spinon this tends to make the hole spectrum more incoherent than the mean fleld 
result. Second as discussed extensively in previous sections the gauge fleld leads to attraction 
between the holons and spinous. This gauge-induced attractive interaction tends to make 
the hole spectrum more coherent than the mean fleld result. These two kinds of effects may 
be viewed as involving self-energy and vertex corrections due to gauge interactions to the 
diagrams that contribute to the electron Greens functions. It is expected that these effects 
do not qualitatively change the mean fleld results to be presented below, at least at the high 
temperatures of interest. 

Turning to the calculation of the spectral function in the cartoon holon-metal phase 
described above, we flnd good comparison with recent temperature dependent photo- 
emission experiments in the normal "pseudo-gap" state of the cuprates^. In the car- 
toon state the propagators for the z and / fleld are, Gi, = —{zz'^) = — 2 and 

Gi^^^ = -{fP) = where is the spinon dispersion and is the kinetic energy of 



Neglecting fluctuations of the gauge fleld and setting A = 0, Eq. (15.31) becomes after 
Fourier-transformation, 



the holons. 




(5.5) 
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FIG. 4: The change in the spectral response as the parameter rj (spectral broadening) is increased. 
In units of the spinon gap T = 0.5. From left to right rj = 0.2, 0.4, 0.8. In the limit of perfect 
resolution the response consists of ellipses centered around the Kp. The spectral peaks are due to 
the higher density of states perpendicular to the zone diagonals due to the larger holon mass in 
this direction. The rest of the results shown in this paper have ry = 0.8 to model the experimental 
resolution. 




FIG. 5: Evolution of A^^(k,0) as a function of T in the first quadrant of the square lattice BZ, 
calculated with Fig. [3]A.. These results bear close resemblance to the shrinking "Fermi Arcs" seen in 
photo-emission experiments on the underdoped cuprates^^. The numerical evaluation of Eq. (j5.6p 
was made with mix/miy = 0.05 to model the anisotropic holon dispersion. These results are in 
the regime where the spectral resolution rj > T, as in current experiments^^. Note that this figure 
has A = and the "Fermi arcs" shrinks to (^, ^) at T = 0. 

The spectral function of the electron is thus a convolution of the spectral functions of its 
composite "holon" and "spinon" particles, diagrammatically this may be simply represented 
as Fig.EKA). 

We can now easily complete the Matsubara sums and analytically continue iun —>■ uj + ir], 
to obtain the expression for the spectral density A = —2'^[G^]. Note that the spectral 
functions do not satisfy the usual sum rule J ^A(uj) = 1 for fermions. This is because we 
are working the in projective subspace of the t-J model. 

Let us consider u = in the "holon-metal" phase, then Eq. (15.50 becomes, 

A^^^ik, 0) = 2vri 5^ ^^^^h^(i?) + n,{e)], (5.6) 

p,v 
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where we have used e = ep-k+Kp and E = Ep for shorthand. Note that at T = this 
evaluates to zero and thus these contributions correspond to propagation of physical electrons 
in the presence of thermally excited spinous and holons. Using = + (for 
simplicity we set the chemical potential of the holons /i = 0) and E^ = a/A^ + c^fc^, we can 
compute the zero energy spectral weight in /c-space. We expect that A{k, 0) should have 
weight only for those momenta, k at which a relative shift by k of the spinon dispersion with 
respect to the holon dispersion causes the bottom of the gapped spinon dispersion to have 
the same energy as the holon. We thus expect spectral weights on ellipses centered around 
the Kp points. We note however that current experiments are carried out in the regime 
where the spectral resolution (16-20meV) is larger than T (110-200 K)^. For this reason, 
the ellipses predicted by our theory can not be resolved, as illustrated in Fig. |H To make 
explicit contact with the experiments, we shall henceforth work in the regime where rj > T 
[f] is the imaginary part used in the analytic continuation). 

The results are displayed in Fig. [5] as a function of the temperature T, they have many 
features similar to the spectral function measured in recent photo-emission experiments^^. 
At finite-T, there are four "Fermi arcs" which have lengths that shrink with lowering T. 
In the holon-metal theory this shrinkage is due to that fact that as T is lowered there are 
fewer and fewer thermally excited spinous present that are required to have an electronic 
excitation. An important difference between the experiments and the results presented in 
Fig. [5] is that in the photo-emission experiments in the pseudo-gap phase, the point to which 
the "Fermi arcs" shrink is shifted away from the Kp. We now show that this observation 
follows on including the effect of Eq. (12.121) . 

We have so far excluded the "Shraiman-Siggia" term that couples the holons and the 
spinous, Eq. fl2.12p . i.e., we have used A = 0. The first order correction to the A'^ due to 
Eq. fl2.12p is given by the diagram shown in Fig. EJ^B). We now focus on the momenta close 
to ii'i = (|, |) and hence drop the three other square lattice symmetry related terms. After 
Fourier transformation we obtain, 



^ 2zA 



2^i{k+ p)xf{p, g, k, iur, 



(5.7) 



where. 



f{p,q,k,iuJn) = g{p,k,iuJn)giq,k,iuJn), 

g{p,k,iujn) = T^Gli^^^Gl+p^i^^+i^^^ (5.8) 
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In order to calculate the spectral function we need to complete the Matsubara sums and then 
continue iun — > u+irj and get the imaginary part of /. Note that because one can always put 
p —>■ —p in the sums above the spectral function contribution changes sign when k^ — > —k^, 
i.e. depending on the sign of A, it shifts the spectral weight in A'^{k, 0) towards(away) from 
the F-point along the zone diagonal. 
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T = 0.5 A, 



T = 0.5 A, 



T = 0.5 A, 





T = 0.1 A, 




T = 0.1 A, 



T = 0.1 A, 




FIG. 6: (color online) The first column shows the shift of the "nodal point" of Fermi arcs from the 
Kp on inclusion of the effect of Eq. (j2.12p at first order in perturbation theory, Fig. [3)3; the plot is 
of A'^ = A'^^^^ +A'^^^\ The next two show the individual contributions, A'^^^^ and A^^^\ The second 
row is with the same physical parameters as the first set but at a lower T = O.IA^. All these data 
are in the regime where r] > T. Lines joining (0, 0) to (vr, vr) and (0, vr) to {it, 0) have been included 
for clarity. 

For = 0, we obtain for the function g{p, k,uj = 0), 



1 

2E 



nB{E) - npie) + I usiE) + npie) 



E + e 



E 



71 



- — [ns{E)+nF{e)] 6{E 



(5.9) 



where we have uses the short-hand, E = E^^p and e = ep. 

We now evaluate the sums in Eq. (15. 7p numerically. Results are presented in Fig. O 
Inclusion of the A term does not change the general features of the shrinking "Fermi arcs" 
observed in Fig. [5l instead it shifts the nodal point away from the Kp, just as is seen in the 
experiments. It should be noted that the holons still have their dispersion minima at the 
Kp, only the physical electrons have their spectral weight shifted away from Kp because of 
A. This should be contrasted with the long-range ordered Neel phase where even Eq. (I2.12p 
cannot change the pinning of the center of the spectral weight away from the Kp (as observed 
in experiments^). 
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VI. CONCLUSIONS 



The starting point of this paper was the deconfined quantum critical point between the 
Neel and VBS ground states of a 5* = 1/2 square lattice antiferromagnet. We added a small 
density of holes to both states and found distinct Fermi liquid states. Their difference was 
characterized in Eq. fll.ll) . and ultimately had its origin in the following results: (i) the 
dispersion of a single hole in a Neel state can have the minimum of its dispersion pinned 
at the symmetrical points Kp (Fig [2]) over a finite range of parameters; (m) the minimum 
of the dispersion of a single hole in the VBS state is generically shifted away from the Kp, 
as established in Eq. (14.91) : the binding of the holon and spinon into a charge e S = 1/2 
quasiparticle in the VBS phase was crucial in obtaining this shift. We also argued that at 
finite hole doping, 6, these Fermi liquid states were likely not connected by a direct quantum 
phase transition, but that the insulating deconfined quantum critical point opened up into 
a non-Fermi-liquid holon metal phase. 

We also discussed the finite temperature, "quantum critical" , holon metal phase obtained 
above the low-temperature ordered phases. The spectra shown in Section IV Al were found 
to be consistent with key characteristics of recent observations of "Fermi-arc" spectra in the 
pseudogap phase of the cuprates. In particular, we showed that the peak in the spectral 
weight was shifted away from the symmetrical Kp position. This shift was linked to the 
same term in the Lagrangian that caused the shift of the dispersion minimum in the VBS 
state, a term originally introduced by Shraiman and Siggia^. 

Finally, we note that our boundary CFT in Section |V] has features that are reminiscent 
of self-consistent impurity models^ of 'local' quantum phase transitions in Kondo lattice 
systems. However, the details are different, and can only be addressed in a quantum field- 
theoretic framework which we have described. This work therefore also serves to place 
studies of quantum transitions in the 'dynamical mean field theory' approacb^i^ in the 
field-theoretic context. 
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APPENDIX A: SINGLE HOLE IN THE U(l) SPIN LIQUID 

In this appendix, we calculate the dispersion of the holons due to absorption and emission 
of the spinous. We will work within the framework of the t — J model Eq. (12.11) . using the 
formalism of Refs. Blssl. 
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FIG. 7: The two diagrams contributing to the holon self-energy in Eq. (|A7p . The sohd hnes are 
the holon propogators and the double-dashed lines are the spinous. 



First consider the second term in Eq. (12. ip . using a Schwinger-boson mean- field theory^i^ 
we arrive at Eq. (12. 5p . After Fourier transformation, this Hamiltonian can be diagonalized 
into the form 

Hj = Y, efc(«L«fc. + PlM (Al) 



ka 



by introducing two new particles, 



aka = Ukbka + Vkh_ka 
Pka = Ukha + Vkb^-ka 



(A2) 



where 



7fc 



J2^^k■6 



COS 



k -\- k 



^ a' ] + cos 



k' - k' 



^ a 



(A3) 



where k'^^y^ G [— ^, ^] (i-e. BZ'), where a' = \f2a. u{v)k are real and u{v)k = u{v)^k- A 
comparison of terms gives, Uk = cosh^^ and = sinh^^, where tanh2^fc = —Q'~fk/^- The 
dispersion of the quasiparticles is = (A^ — (Qta:)^)^^^- Using these we can calculate 
all the propagators of interest. 



— {bkab-ka) 



■-— -2 =^?(/c,^Z/„) 

Qlk 



= T{k, iUri 



(A4) 
(A5) 



We now define a = Q/X and use this parameter instead of Q. Then, e.g., the dispersion 
becomes = A(l — (0:7^)^)^/^, where < a < 0.25. 

The hopping term in the t — J model in terms of Schwinger bosons after Fourier trans- 
formation gives, 

Ht = t ^ 'jkS era' 

fij (A6) 



Working to second order in the hopping, we now calculate the self-energy of the /± 
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particles due to the Ht- There are two diagrams at order (shown in Fig. [7]) that need to 
be summed up, they evaluate to, 

T.^{k,iUn) = [2]t^^^-fg_pg{iUra,q)G{iuJn-iuJi + iiym,k + q-p)Qf{iuJi,p) (A7) 
T,^{k, iUn) = [-2]t'^^ ^ -fq_p'yg+k^{iurn, q)^{i^n " i^i + ii^m, k + q- p)gf {iuji,p) 

The factors in [...] appear from the spin sums. The sum of the above diagrams becomes (we 
assume here that A{p,uj < 0) = 0), 

S^(it^„, k) = '^[t'^{-fg.pUk+q-pVg - Jq+kUqVk+q-pf]gf {iuJn " " Cfe+g-p, k) (A8) 

g,p 

A full self-consistent solution of the above equation corresponds to summing all diagrams 
that have non-crossing boson lines. The dispersion of the f-particles is obtained by looking 
for the poles of the interacting Greens' function, i.e. Uk = ^[T,{uk, k)]. We shall not 
undertake the fully self-consistent calculation here, instead we simply use the self energy to 
second-order in t which gives. 

q,p ^k+q-p eg 

where Ft2(^k,p,q) = t'^^'jq-pUk+q-pVq — '~^q+kUqVk+q-pY ■ This result is plotted in Fig. [HI As 
discussed in the text, the holon disperion has minima at the (ib7r/2, ib7r/2). The band masses 
at the minima are highly anisotropic with the holon dispersion much lighter along the zone 
diagonal than perpendicular to it. 



APPENDIX B: CRITICAL PROPERTIES OF THE CP^^ MODEL 

As argued elsewhere, and reviewed in Section [H the critical theory of the Neel-VBS 
transition on the square lattice is described by the CP^^^ model. Here we will review 
the expansion of this theory, obtain an effective field theory suitable for bound state 
problems, and compute a scaling dimension in the boundary CFT associated with a charged 
impurity. 

We are interested in the properties of the field theory 
Z = j VzaV\VAf,exp{-S) 

(Bl) 



S = I d\d 



^{\{d^-lA^)Zo.\^ + zX{\Za\^-N)} 

Here the index a = 1 . . . A^, and r is imaginary time. The theory has a conformally invariant 
critical point at a = g^, and we want to describe the spectrum for g > g^.. 
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FIG. 8: (color online) The dispersion ojkt^/X obtained from Eq. ()A9p . The points Kp where the 
holon dispersion have their minima are marked for clarity. We have used a = 0.24 here, note 
however that these results are generic. 

First, let us review the properties at = oo. The quantum critical point is at 5^ = Qc-, 
where 

''''' ' (B2) 



9c Jo 87r3 p2 

Here, A is an ultraviolet cutoff, and we will only be interested in phenomena at scales much 
smaller than A. 

For g > Qc, we have a large N solution for the spin- gap U(l) spin liquid with (za) = 0, 
iX = Aq, where 

Here Aq is a measure of the spin gap at = 00; we will discuss the renormalization of this 
to the energy scale A in the following subsection. 

For g < gc, we have the Neel phase, which breaks the SU(A^) rotation symmetry with 
{za) = mo\/N6a,i and iX = 0. The value of mg is given by 

The properties of this phase are easy to understand by spin-wave theory ( "chiral perturbation 
theory"), and we won't consider it further. 

We now consider the structure of fluctuations for g > gc. The 1/N expansion can be 
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easily generated by a Feynman graph expansion of the theory 
Z ^ J VzaVXVAi,exp{-S) 

S = J (frdT[{\{d^-lA^)Za\^ + Al\Za\^+lX\Za\^}] 



(B5) 



+ 



8^ 



Here we have rescaled z^ by ^/g, and there is no further dependence upon g. The kernels n(p) 
and K{p) represent one-loop z^ contributions. We avoid double counting in the Feynman 
graph expansion simply by avoiding any further bubble self-energy contributions to the A 
and propagators. The explicit values of the kernels are 



n(p) 



87r3 (g2 + A2)((p + g)2 + A, 



1 ,-1 /2Ao 
cot ' 



Attp \ p 

l/(87rAo) p < Ao 
l/(8p) p > Ao 



(B6) 



and 



Kip) 



-/ 

i: 



d^q 



<f-{p- <lf/P^ 



87r3((g+p/2)2 + A2)((g-p/2)2 + Ag) 



+ 2 



dx 
2^ 



p^xil - x) + A§ - Ao 



87rp 

pV(247rAo) 
p/16 



2Ao 



P 

P < Ao 
P>Ao ' 



Ao 
47r 



d^q 



87r3 g2 + Ag 



(B7) 



where we have used the Feynman parameter method to evaluate the integrals. 

The following subsection will formulate an approach suitable for bound state problems, 
and examine the response to an impurity. 



1. Non-relativistic effective field theory 

We are interested here primarily in the physics at scales p < Aq. As in ordinary quantum 
electrodynamics, this regime is conveniently expressed in terms of a non-relativistic effective 
field theory. Wc will sec that the low-velocity expansion associated with the effective field 
theory is justified because the typical velocity in a bound state scales as 1/y/N. 
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We first need to express the relativistic theory in (IBSP in terms of the non-relativistic 
particle and anti-particle eigenstates. To this end, we decouple the kinetic term for the Za 
by a Hubbard-Stratanovich field IIq,: 

\{dr - iA^)z^\^ |n„p + ilil{dr - iA^)za + c.c. (B8) 



Then we parameterize 



^ {Pa + hi) 



n. = {pa - hi) (B9) 



where pa and /i" will be the spinon and anti-spinon excitations. For the SU(2) case, with 
= 2, we will make the replacement /i" e'^^hp, so that both the spinons and anti-spinons 

transform under the same representation of SU(2). 

We now insert flB9P in the effective action flB5l) . and expand in powers of spatial gradients. 

The structure of this expansion determines the form of the effective field theory. With some 

additional approximations which will shortly be justified, we postulate, up to order 

the following non-relativistic effective field theory 



Znr = J VpaVhaVArexpi-SNR 
Snr = I (frdr 



(BIO) 



Here we have defined the renormalized energy scale A so that the co-efficient of the elec- 
trostatic term {WAtY term is A^/(487rA). At = 00, we have A = Aq; there will be 
corrections at order 1/A^, which can be computed from fIBSp . which will depend upon the 
cutoff A, and which ensure that A ~ ((7 — g^)^ ^ with4^ z/ = 1 — 48/(7r^A^). The numerical 
constant Ci is a universal number to be determined by a matching computation to the fully 
relativistic theory (IBSD : we will not carry out such a computation here. 

The dominant interaction is the 'Coulomb' interaction between the pa and the /i", me- 
diated by the electrostatic potential Ar. This leads to a potential energy between opposite 
unit charges of the form 

247rA f d^k{l-(^^'' 

V[r) = 



N J 47r2 A;2 
12A 

= ^[ln(rA) + C2] . (Bll) 
Here (2 appears to be a second unknown constant, but we have the freedom to set C2 = 
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without modifying any observable property. This is because Z^v^ apphes only to spinon/anti- 
spinon configurations which have net charge zero, and for these C2 can be absorbed into a 
redefinition of (i. 

From this potential we see that the spinon and anti-spinon form a bound state with 
a spatial extent ~ a/ZV/A. This means that the velocity expansion holds in powers of 
1/\/N. One can also see from this that the coupling to the spatial components of the 
gauge field, A, lead to corrections which are higher order in from the terms included. 
Short-range interaction terms (induced from A fiuctuations) like are also easily 

seen in perturbation theory to be higher order: they induce corrections which depend upon 
the value of the bound state function at the origin ~ 1/^/N. These arguments justify our 
dropping A and A fiuctuations in Znr. 

2. Boundary exponent of an impurity 

As motivated in Section |Vl we consider here the problem of a static charged impurity 
coupled to the 2+1 dimensional critical point of CP^~^ model. For this we modify the 
theory in Eq. fIBip to 



Z = j VzJ)fVXVAf, exp (^-5^ 
S = j d^drC, + J drp + ^0 - iA^ir = 0, r)^ /, (B12) 



l^z = \{d^-ieoA^)Za\'^ + s\Za\^ + — {\Za\'^y 



Note that the Grassman field / depends only upon r, while the bosonic fields vary over both 
r and r. For completeness, we add "bare" terms corresponding to the spinon interaction and 
the gauge field kinetic energy, and we will see that the critical exponents are independent 
of these terms. We will now compute the scaling dimension r]h/2 of the composite operator 
h{T) = z^ir = 0,T)p{T). 

At the critical point s = Sc, the bare boson propagator is Go{k) = The dressed 

photon propagator can be obtained by summing over ^-bubbles and has the form 

D,u{p) = + IP)-' (5,. - C^) , (B13) 

where ( is the arbitrary "gauge parameter" and 

7 = Nel/16. 
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In a similar way, one can obtain the dressed 2;-field interaction vertex 



^ ' l + 2uoNU{ky 

where n(/c) = l/(8/c) is the bosonic bubble. The large- expansion is thus the expansion 
in dressed photon lines and dressed ^-vertices. 

The anomalous dimensions 77^ and ?7/ of the fields z and / can be read off as the coefficients 
in the contributions of the type —rjzp^lYiij)) and irjfUj\Yi{uj) to the self-energies ^^(p) and 
S/(u;) which enter the corresponding dressed propagators G~^{p) = p'^ + T.z{p) and G'J^(cj) = 
—iuj + ^f{uj) (here u is counted from the reference energy Eq). To the order of the 
self-energies are given by 

- el{2p^ - k^){2p, - k,)D^,{k)^, (B14) 

/Sk 1 



and after tedious but straightforward calculation one obtains 

^^ = ^-#^-(1-0;.^, (B16) 

% = -(l + C)^, (B17) 

One should realize that, taken separately, rj^ and r^/ are both gauge-dependent and do not 
have any direct physical sense. It turns out that there is no contribution to rjh from the 
vertex diagram describing an exchange of a photon between the z and / particles, so to the 
1/N order, the scaling dimension rih/2 of h is given by 1/2 + rj^ + rjf, which in turn yields 
the (gauge-independent) expression (15. 2p for the boundary exponent rjh. 
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